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A method to determine the confidence limits for the abscissa of the intersection of two linear 
regressions has been developed. This method does not require the assumption of equal variance for 
the two regressions, as was necessary with previous methods. A numerical example is included on 
thermodynamic, glass transition data for which this method is applicable. Comparisons are made 
between the results using equal and unequal variance assumptions. A FORTRAN subroutine is included 
for computations using both assumptions. 
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1 . Introduction 

In this paper confidence limits are derived for the abscissa of the intersection point of two 
linear regressions. As illustrated in figure 1, it is assumed that there are independent realizations 
of two linear regressions at ordered values of a single independent variable X. The problem is to 
determine confidence limits for the abscissa Xo of the intersection point of the two linear regressions. 

This problem was originally considered and solved by Fisher [l] 1 and later by Kastenbaum 
[2], but under the assumption of a common error variance for the two regression lines. Hinkley [3] 
derived the asymptotic distribution of the maximum likelihood estimator for the intersection point 
abscissa of two linear regressions, also assuming equal variances. Robison [4] presented estimators 
and confidence limits for the intersection point abscissa of two polynomial regressions under 
the equal variance assumption. Hudson [5] derived least squares estimators (but no confidence 
limits) for the intersection point abscissas of two or more general regression models for the unequal 
variance case. 

The principal contribution of the present paper is a derivation of confidence limits for the 
intersection point abscissa without resorting to the equal variance assumption. Although only 
the linear regression case is considered below, extension of results to more complicated regression 
models is possible. 

The results of this derivation have direct application to thermodynamic data to evaluate the 
uncertainties in the second order transition temperature T 2 and the glass transition temperature T g . 



AMS Subject Classification: 62.55. 

1 Figures in brackets indicate the literature references at the end of this paper. 
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FIGURE 1. Intersection of two linear regressions. 



The second order transition is an equilibrium phenomenon, whereas the glass transition is a non- 
equilibrium one for which T 9 depends upon the thermodynamic history by which the glass is formed. 
Both transitions may be manifested by discontinuities in the second order properties (for example, 
thermal expansivity or heat capacity) with respect to temperature. The corresponding first order 
properties (volume or enthalpy) vary with temperature in a manner similar to Y versus X in figure 1 
for which the transition temperature may be determined by the intersection of the two straight lines 
representative of the response in the two regions. Accordingly, this method may be applied directly 
to any of the first order properties to determine the confidence bands of the transition temperature 
as defined 2 from the intersection of the two regression lines. The method is sufficiently general that 
it is applicable to cases where data are not available in the proximity of the intersection, which, 
of course, leads to a less accurate prediction of its abscissa. 

Two methods are described to evaluate the confidence limits for the abscissa of intersection. 
These are based on the assumptions of equal and unequal variances for the regression lines deter- 
mined from the two sets. The advantage of the equal variance case is that the underlying assumption 
of equal variance permits a relatively simple development of the method, whereas the advantage of 
the unequal variance case is, of course, in its more general application. The ANSI FORTRAN sub- 
routine to calculate the confidence limits for both cases is given in the appendix, and numerical 
examples are included in later discussions. 



2. Derivation of Confidence Limits 

A A 

As illustrated in figure 1, let Yi(X) = a x + b x X and Y 2 (X) = a 2 + b^ be two fitted least 
squares regressions based on ri\ and n 2 observations (Xij, i= 1, 2;j= 1,2, . . ., m) respectively, 
where the random variables Yi(X) are assumed to be normally and independently distributed about 
the true unknown regression lines ^i{X) = ai~\- fitX with variances cr?. (ai, fit, and a* are true, 
but unknown, parameter values.) 

The problem is to determine confidence limits for the abscissa Xq of the intersection point 
(X , Y ) of the two unknown regression lines. 

The starting point for our derivation is the fact that a,- and 6*, the usual least squares esti- 
mators, are linear combinations of normally distributed variables and hence themselves are 
normally distributed. Thus for all X, Yi(X) —^cn + b[X is also normally distributed, as is the 
difference Y\(X) — Y 2 (X). Since Yi(X) — Y 2 (X) has mean 1*\{X) — /J, 2 (X) and variance 



- This definition of the glass temperature is one of several which are not necessarily equivalent. 
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a 2 t Vt(X)+crlV 2 (X), where 



V i (X) = ±+ {X J i} \ £=1,2 



1 n i 
^ = ~2^ U ' '"1.2 (1) 






then it follows that all X 



[YdX)-Y 2 (X)]-[MX)-^(X)] 
VcrfF,(J) + <rlV 2 (X) 



(2) 



has a standardized (mean 0, variance 1) normal distribution. 

Consider first of all the equal variance case (a 2 = cr\ = cr 2 , say), treated earlier fl]. It is seen 
that eq (2) simplifies to 

[Y^X) - 9,00] - \jii(X) - »,(*)] 

<tVV,(X)+ V 2 (X) (3) 

Thus 

[Y^-YAX)] - [MX) ~ H*(X)] 



s VV,(X) + V 2 (X) 



(4) 



has a t distribution with ti\ + ra 2 — 4 degrees of freedom, where 5 is the square root of the pooled 
unbiased estimator of a 2 : 



/ (n, -2)s» + (m 
' n, + n> — 



n >- 2 ^, (5) 



and where s 2 is the unbiased estimator for cr 2 given by the residual variance obtained from the 
least squares fit of line i. 

The quantity (4) is distributed as t for all X, and in particular for X = Xo in which case 
IXi(X) — fJLziX) vanishes. It then follows that for every probability p(0 < p < 1) 

p,.b(( ftt*>-ftc»»> )'«( c /i±£))n- p (6 , 

where G t f — - — J is the 100 ( — - — j percent point of the t distribution with n t + n 2 — 4 degrees of 

freedom. The squared numerator [Yi(Xo)— Y 2 (Xo)] 2 and squared denominator s 2 [Vi(Xo) 
-f V 2 (Xo) ] above are both quadratics in X . Thus eq (6) is of the form 

Prob I jr- ^ G 2 1 (or equivalently Prob {Qi — G 2 Q 2 ^ 0}) where Qi denote quadratic expressions, 
and G denotes G t I — - — j. Since Q\ — G 2 Q 2 is itself a quadratic (say AX 2 + BXo + C), then the 
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desired confidence interval consists of those X such that AX 2 Q + BXo + C « and the desired con- 
fidence limits for Xo are simply the two roots of AX* + BXo + C = 0. 

The coefficients A, B, and C are rather tediously derived; omitting the intermediate algebra, 
these coefficients are 



A= (&,-o 2 ) 2 -s 2 G 2 



L s 2 s 2 



B = 2(a, - a 2 ) (6, - b t ) + 2s 2 G 2 



C= (a,-a 2 ) 2 -s 2 G 2 



± + ± + xi , xr 



? + s 2 J 
|| ||] 

S 2 S 2 J' 



(7) 



It is of interest to note that the midpoint — B/2A of the confidence interval is not identical to the 
maximum likelihood estimator 



X = - 



a\ — a,2 



(8) 



however, the two estimators converge to a common limit as s 2 — » 0. For a less detailed derivation of 
the above confidence interval results, see references [1-3]. 

For the more general cr\ ^ cr\ case, the distributional complications which arise are identical 
to those encountered in the Behrens-Fisher 2-means problem [6]. The solution outlined below is 
an extension of the Welch-Aspin [7, 8, 9] solution to the Behrens-Fisher problem. Proceeding in a 
manner analogous to Welch and Aspin, we see that for allZ, 



[Y l (X)-Y 2 (X)]-[ IJ . i (X)- l j, 2 (X)] 



(9) 



has approximately a t distribution with v degrees of freedom, where v is given by 



1 



Vi(X ) 



1 



v n 1 -2\V l (X )+V 2 (X, 



o)/ n 2 -2\V 1 (Xo) + V 2 (X (> )) 



(10) 



A minor complication at this point is the fact this Xo in (10) is an unknown parameter. In order 
to arrive at an approximate value for v the maximum likelihood estimate Xo (as defined in (8)) is 
used. The development now parallels the a 2 = al case, and so again 100 p confidence limits for 
Xo are obtained by solving a quadratic equation, but with the following slightly-modified coefficients: 



^ (4l _ w ,_c[£fj|] 



B = 2(a 1 - a2 )(6 1 -W-r2G 2 [^i^] 



(11) 



C=(a 1 -a 2 ) 2 -G 2 



Tii n 2 S? 



J? si J 
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Note that although the unequal variance case basically parallels the equal variance case, the 
former does not numerically reduce to the latter when sj = s\. The confidence intervals usually 
are distinct even when s 2 { = s§. In deciding whether to use the equal variance results or the 
unequal variance results, there are two considerations to be taken into account: physical theory 
and the variance ratio F test. The general rule of thumb is to always use the unequal variance 
results unless the following two conditions are simultaneously satisfied: (1) there is a firm physical 
reason for assuming equal variances; and (2) the F ratio Sjf/Sf is not significantly different from unity. 
This recommended procedure is at times conservative, but more important, it is less assumption- 
dependent and hence more generally valid. 

In the appendix, an ANSI FORTRAN subroutine is presented which first of all tests to see 
whether the equal variance assumption is tenable, and secondly computes confidence limits for 
Xo for both the equal variance and the unequal variance cases. The subroutine is double precision 
in input, internal operation, and output. Definitions of input parameters and instructions for the 
use of the subroutine are given in the comment statements within the subroutine. 

3. Example: Calculation of the 95 percent Confidence Limits for the 
Pressure-Dependent Glass Temperature 

As an example, the 95 percent confidence limits have been calculated for the pressure- 
dependent glass transition temperature of poly( vinyl acetate), an amorphous, glass-forming polymer. 
The PVT (pressure- volume-temperature) data for the glass, designated as set 1, are given in 
table 1, and the data for the liquid, designated as set 2, are given in table 2. Both sets of data, 









Table 1. 


Specific volu 


me values for 


the glassy regi 


Ion 






T/P 





100 


200 


300 


400 


500 


600 


700 


800 


-30 


0.82884 


0.82654 


0.82438 


0.82225 


0.82017 


0.81816 


0.81613 


0.81419 


0.81237 


-25 


.83002 


.82768 


.82547 


.82334 


.82121 


.81914 


.81715 


.81516 


.81327 


-20 


.83110 


.82874 


.82648 


.82430 


.82218 


.82010 


.81808 


.81606 


.81415 


-15 


.83238 


.82997 


.82767 


.82549 


.82331 


.82117 


.81919 


.81717 


.81522 


-10 


.83350 


.83117 


.82886 


.82663 


.82443 


.82224 


.82022 


.81817 


.81620 


-5 


.83467 


.83233 


.83003 


.82774 


.82557 


.82340 


.82131 


.81933 


.81729 





.83581 


.83351 


.83115 


.82889 


.82666 


.82447 


.82237 


.82033 


.81825 


5 


.83706 


.a3478 


.83241 


.83012 


.82794 


.82569 


.82355 


.82140 


.81941 


10 


.83812 


.83593 


.83359 


.83130 


.82903 


.82679 


.82465 


.82248 


.82040 


15 


.83937 


















20 


.84061 



















Units: T in °C, P in Bars, v in cm 3 /g. 



Table 2. Specific volume data for the liquid region 



TIP 





100 


200 


300 


400 


500 


600 


700 


800 


35 


0.84572 


0.84148 
















40 


.84870 


.84444 


0.84041 


0.83670 












45 


.85174 


.84733 


.84321 


.83921 


0.83549 


0.83184 








50 


.85486 


.85041 


.84608 


.84206 


.83817 


.83447 


0.83087 


0.82757 




55 


.85791 


.85349 


.84894 


.84491 


.84092 


.83709 


.83340 


.82988 


0.82666 


60 


.86104 


.85628 


.85179 


.84769 


.84367 


.83980 


.83607 


.83256 


.82904 


65 


.86407 


.85933 


.85472 


.85052 


.84648 


.84250 


.83874 


.83523 


.83174 


70 


.86723 


.86218 


.85755 


.85324 


.84913 


.84507 


.84108 


.83734 


.83378 


75 


.87038 


.86536 


.86038 


.85594 


.85171 


.84762 


.84367 


.84000 


.83626 


80 


.87343 


.86829 


.86342 


.85881 


.85453 


.85036 


.84641 


.84257 


.83889 


85 


.87669 


.87140 


.86636 


.86169 


.85728 


.85308 


.84911 


.84532 


.84155 


90 


.87986 


.87438 


.86923 


.86444 


.86000 


.85574 


.85172 


.84776 


.84396 


95 


.88301 


.87741 


.87207 


.86722 


.86269 


.85823 


.85424 


.85025 


.84641 


100 


.88622 


.88043 


.87500 


.87012 


.86534 


.86095 


.85678 


.85284 


.84904 



Units: T in °C, P in bars, v in cm'Vg 
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taken from reference [10J, are illustrated in figure 2, where the specific volume v taken from the 
above tables is plotted as a function of temperature at the different pressures indicated. The 
specific volume data Vi(T, P) for the glass (at temperatures between —30 and 20 °C) were obtained 
after the glass was formed at atmospheric pressure at a constant cooling rate of 5 °C/h commencing 
at equilibrium in the liquid region. The data Vi{T, P) for the liquid (35 to 100 °C) are taken to 
be in true equilibrium. A characteristic of a glass is its large increase in its viscoelastic relaxation 
times with decreasing temperature. The glass transition is a manifestation of the increase of 
relaxation times during formation of the glass. Although the data shown here for the glass are 
not representative of the attainment of true equilibrium, the relaxation times at these temperatures 
are very long in comparison to experimental times. Accordingly, the glass data may be treated as 
being in apparent equilibrium (pseudoequilibrium) over any practical time scale, as has been 
done in this case. 

Since the volume isobars vary nearly linearly with temperature in both the liquid and glass 
regions, and since the slopes for each isobar in these regions are distinct, the results of section 
2 are applicable. By treating each of the two regions separately, each of the volume isobars may 
be represented by two distinct regression lines of the form 



ii = d\ + 6,T, i = 1, 2. 



(12) 



The intersection of these two regression lines at each pressure gives the ''maximum likelihood 
estimator" 



r„=- 



at— q 2 
6]— b 2 



(13) 



which is taken here as an estimate of the true, but unknown, glass temperature 7^for the particular 
thermodynamic history by which the glass was formed. In figure 2 the regression lines of the 
form of eq (12) are shown by the solid lines with positive slope at each pressure. The solid line 




-30 -20 -10 



90 100 



FIGURE 2. Specific volume of polyvinyl acetate) as a function of temperature at 

different pressures. 

An estimate of the glass transition is obtained from the intersection of the volume isobars taken to be linear with 
temperature. 
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with negative slope gives the positions of T g , which depend strongly on the pressure. It is seen 
that some of the data representative of the liquid region fall to the left of T u . Although this result 
may b^ paradoxical, it is simply a consequence of the fact that the liquid data are equilibrium 
data, and the glass are pseudoequilibrium data 4 resulting from the particular thermodynamic 
history by which the glass was formed. 

Ninety-five percent confidence limits for T g were calculated using both the equal and unequal 
variance assumptions. The paired statistical quantities necessary to calculate the confidence limits 
are given in tables 3 and 4 for the glass and liquid, respectively, where P is the pressure, a\ is the 
intercept, bi is the slope, st is the standard deviation of the residuals about the regression lines, 

TABLE 3. Statistical data for the (Glassy Region 



p 


a x 


6, X 10 3 


5l X 10 4 


n x 


f, 


s\ 





0.835843 


0.23469 


0.495 


11 


-5.0 


2750.0 


100 


.833531 


.23590 


.565 


9 


-10.0 


1500.0 


200 


.831192 


.23123 


.778 


9 


-10.0 


1500.0 


300 


.828928 


.22660 


.855 


9 


-10.0 


1500.0 


400 


.826716 


.22270 


.910 


9 


-10.0 


1500.0 


500 


.821512 


.21713 


1.021 


9 


-10.0 


1500.0 


600 


.822117 


.21327 


.793 


9 


-10.0 


1500.0 


700 


.820328 


.20860 


.737 


9 


-10.0 


1500.0 


800 


.818301 


.20270 


.897 


9 


-10.0 


1500.0 



TABLE 4. Statistical data for the Liquid Region 



p 


a 2 


6 2 X 10 :J 


5 2 X 10 4 


n 2 


fi 


si 





0.823686 


0.62355 


1.146 


14 


67.5 


5687.5 


100 


.820417 


.59934 


.869 


14 


67.5 


5687.5 


200 


.817195 


.57760 


.683 


13 


70.0 


4550.0 


300 


.814189 


.55834 


.832 


13 


70.0 


4550.0 


400 


.811039 


.54362 


.565 


12 


72.5 


3575.0 


500 


.808010 


.52946 


.589 


12 


72.5 


3575.0 


600 


.804836 


.51987 


.952 


11 


75.0 


2750.0 


700 


.802096 


.50698 


1.341 


11 


75.0 


2750.0 


800 


.799213 


.49703 


1.564 


10 


77.5 


2062.5 



in is the number of observations, T\ is the average temperature, and S? is the sum of the squared 
temperature deviations: 

S^^iTij-Ti)*, j=i,2. 

The results of these calculations are summarized in table 5 for both the equal and unequal 
variance cases. G is the 97.5 percent point of the t v distribution with v degrees of freedom where 
v = rii + riz— 4 in the equal variance case, and where v is defined by eq (10) in the unequal variance 
case. —A/2B is the midpoint of the confidence interval where A and B are defined by eqs (7) in 
the equal variance case and eqs (11) in the unequal variance case. Aris the 95 percent confidence 
interval width which is equal to the difference of the two roots of the quadratic equation for which 
the constants are defined in eqs (7) and (11) for equal variance and unequal variance, respectively. 
s 1/5-2 is the ratio of residual standard deviations, and R is the ratio of AT" for the equal variance 
case to A7 1 for the unequal variance case. Perusal of table 5 reveals several interesting points. 
First, it is noted that the corresponding midpoints —B/2A are essentially identical for the equal 



4 In principle, at least at temperatures over a limited range below T„ shown here, the specific volume values for the glass will relax to extrapolations of the liquid 
isobars shown provided sufficient time is allowed. Since this time may be as long as several million years (depending on the temperature and pressure), the attain- 
ment of true equilibrium in all cases (as with the liquid values shown) is not practical. 
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Table 5. Estimated values ofT n and 95 percent Confidence Intervals 









Equal variance 






Unequal 


variance 








p 


T 9 


















5 2 


R 






















V 


G 


-B/2A 


AT 


V 


G 


-B/2A 


AT 









31.26 


21 


2.08 


31.27 


0.92 


17 


2.11 


31.26 


0.79 


0.432 


1.162 


]00 


36.08 


19 


2.09 


36.09 


1.18 


9 


2.26 


36.09 


1.02 


0.650 


1.158 


200 


40.41 


18 


2.10 


40.42 


1.26 


9 


2.26 


40.42 


1.44 


1.138 


0.874 


300 


44.43 


18 


2.10 


44.44 


1.62 


8 


2.31 


44.44 


1.80 


1.028 


0.898 


400 


48.85 


17 


2.11 


48.86 


1.56 


8 


2.31 


48.87 


2.07 


1.609 


0.752 


500 


52.84 


17 


2.11 


52.85 


1.84 


8 


2.31 


52.86 


2.53 


1.734 


0.729 


600 


57.34 


16 


2.12 


57.36 


2.24 


8 


2.31 


57.35 


2.21 


.833 


1.014 


700 


61.10 


16 


2.12 


61.13 


3.03 


7 


2.36 


61.11 


2.34 


.549 


1.296 


800 


64.85 


15 


2.13 


64.89 


3.77 


7 


2.36 


64.87 


3.02 


.573 


1.248 



and unequal variance cases, and these, in turn, are essentially identical to the maximum likelihood 
estimator T g . As mentioned earlier these quantities converge to the common limit T g as sf — > 0. 
Secondly, the confidence limits vary considerably between the equal and unequal variance cases 
(the ratio R being at times both greater than or less than unity). In the case where si and 52 are 
nearly equal (P = 300 bars), R differs considerably from unity; and in the case where R is nearly 
unity (P = 700 bars), si and 52 are considerably distinct. This behavior illustrates the fact that 
the equal variance case is not a special case of the unequal variance case as might be expected. 
Thirdly, there appears to be little correlation between the ratio of confidence interval widths and 
the ratio of residual standard deviations. 

The values of the maximum likelihood estimator T g and those for the 95 percent confidence 
interval width AT (each centered about its midpoint — B/2A) for both the equal variance (cr^ = cr\) 
and the unequal variance {cr\ # &\) cases are shown on figure 3 as functions of pressure. The 
median solid line is a quadratic regression of T g with respect to pressure. The two extreme solid 
lines are quadratic regressions of the confidence limits with respect to pressure for the unequal 
variance case only. It is clear that a trend exists between the confidence interval and pressure. 
The principal reason for the divergence of the confidence interval with increasing pressure is 
that the difference (T g — 7\) between the estimated glass temperature T g and the average tem- 
perature for the glass 7\ also increases considerably with pressure as seen from figure 2. Thus 
is apparent the importance of obtaining data in the vicinity of the glass temperature (or, in general, 
the intersection) in order to obtain an accurate estimate. 




Figure 3. Maximum Likelihood Estimator, T g , and 95 percent Confidence Intervals as 

functions of pressure. 

The solid lines are quadratic regressions of these quantities. 
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4. Appendix 
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SUBROUTINE ABSCIS ( AN1 , AN2 , XBAR1 , XBAR2 , SSQDX1 , SSQDX2 , ALP1 , ALP2 , BET1 
1,BET2,RSD1,RSD2) 

THE PURPOSE OF THIS ROUTINE IS TO COMPUTE 95 PERCENT CONFIDENCE LIMITS 

FOR THE ABSCISSA OF THE INTERSECTION POINT OF TWO LINEAR REGRESSIONS 

FOR BOTH THE EQUAL VARIANCES CASE AND THE UNEQUAL VARIANCES CASE 

AN1 = THE NUMBER OF (X,Y) PAIRS FOR LINE 1 

AN2 = THE NUMBER OF (X,Y) PAIRS FOR LINE 2 

XBAR1 = THE SAMPLE MEAN OF THE X VALUES FOR LINE 1 

XBAR2 = THE SAMPLE MEAN OF THE X VALUES FOR LINE 2 

SSQDX1 = THE SUM OF THE SQUARED DEVIATIONS OF THE X VALUES FOR LINE 1 
ABOUT XBAR1 

SSQDX2 = THE SUM OF THE SQUARED DEVIATIONS OF THE X VALUES FOR LINE 2 
ABOUT XBAR2 

ALP1 = THE Y INTERCEPT OF LINE 1 

ALP2 = THE Y INTERCEPT OF LINE 2 

BET1 = THE SLOPE OF LINE 1 

BET2 = THE SLOPE OF LINE 2 

RSDI = THE RESIDUAL STANDARD DEVIATION OF LINE 1 

RSD2 = THE RESIDUAL STANDARD DEVIATION OF LINE 2 

THIS ROUTINE IS DOUBLE PRECISION IN INPUT, INTERNAL OPERATION, AND OUTPUT 

DIMENSION Z(30) 

DOUBLE PRECISION AN1 , AN2 ,XBAR1 ,XBAR2 , SSQDX1 ,SSQDX2 

DOUBLE PRECISION ALP1 , ALP2 ,BET1 ,BET2 ,RSD1 ,RSD2 

DOUBLE PRECISION Z 

DOUBLE PRECISION P,S ,DELALP ,DELBET , AML 

DOUBLE PRECISION ANUEQ , ANUNEQ, GEQ ,GNEQ, PP ,HPP 

DOUBLE PRECISION RATIO , AM , AN, HOLD 

DOUBLE PRECISON Al , A2 , A3 , A, Bl ,B2 ,B3 ,B , CI , C2 , C3 , C 

DOUBLE PRECISION DISC ,ROOTl ,R00T2 , AMID ,WIDTH1 ,WIDTH2 ,RATWID 

DOUBLE PRECISION DSQRT 

DATA Z(1),Z(2),Z(3),Z(4),Z(5),Z(6),Z(7),Z(8),Z(9),Z(10),Z(11), 
1Z(12) ,Z(13) ,Z(14) ,Z(15) ,Z(16) ,Z(17) ,Z(18) ,Z(19) ,Z(20) ,Z(21) , 
1Z(22) ,Z(23) ,Z(24) ,Z(25) ,Z(26) ,Z(27) ,Z(28) ,Z(29) ,Z(30 )/12 .7026D0 , 
14 . 3027D0 , 3 . 1824D0 , 2 . 7764D0 , 2 . 5706D0 , 12 . 4469DO , 2 . 3646D0 , 2 . 3060D0 , 
12 . 2622D0 , 2 . 2281D0 , 2 . 2010D0 , 2 . 1788D0 , 12 . 1604DO , 2 . 1448D0 , 2 . 1315D0 , 
12 . 1199D0 , 2 . 1096D0 , 2 . 1009D0 , 2 . 0930D0 , 12 . 0860D0 , 2 . 0796D0 , 2 . 0739D0 , 
12 . 0687D0 , 2 . 0639DO , 2 . 0595D0 , 2 . 0555DO , 12 . 0518DO , 2 . 0484DO , 2 . 0452D0 , 
12.0423D0/ 



*RSD2*RSD2 ) /( AN1 + AN2-4 . 



COMPUTE NEEDED CONSTANTS 

P=.95D0 

S=DSQRT ( ( ( AN1-2 . ODO ) *RSD1*RSD1+ ( AN2-2 . ODO 1 

1DO)) 

DELALP=ALP1-ALP2 

DELBET=BET1-BET2 

COMPUTE THE MAXIMUM LIKELIHOOD ESTIMATOR 

AML=-DELALP/DELBET 

COMPUTE THE APPROPRIATE DEGREES OF FREEDOM FOR THE T PERCENT POINTS 
FOR THE EQUAL VARIANCES CASE AND FOR THE UNEQUAL VARIANCES CASE 

ANUEQ=ANl+AN2-4 . ODO 

Vl= ( 1 . ODO/AN1 ) + ( AML-XBAR1 ) * ( AML-XBAR1 ) /SSQDX1 
V2= ( 1 . 0D0/AN2 ) + ( AML-XBAR2 ) * ( AML-XBAR2 ) /SSQDX2 

ANUNEQ=(1.0D0/(AN1-2.0D0) ) *( (V1/(V1+V2) ) **2 )+( 1 . ODO/( AN2-2 . ODO ) )*( 
1(V2/(V1+V2))**2) 
ANUNEQ=1 . ODO/ANUNEQ 
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NUEQ=ANUEQ+0 . 1D0 

NUNEQ=ANUNEQ+0 . 1DO 

IF ( NUEQ . LE . 30 ) GEQ=Z ( NUEQ ) 

IF ( NUEQ . GT . 30 ) GEQ=1 . 96D0*DSQRT ( ANUEQ/ ( ANUEQ-2 . ODO ) ) 

IF ( NUNEQ . LE . 30 ) GNEQ=Z ( NUNEQ ) 

IF ( NUNEQ . GT . 30 ) GNEQ=1 . 96D0*DSQRT ( ANUNEQ/ ( ANUNEQ-2 . ODO ) ) 

PP=(P+1.0D0)/2.0D0 

HPP=100.0D0*PP 
C 

C WRITE OUT THE INPUT DATA 

C IT IS ASSUMED THAT THE OUTPUT UNIT IS 6 . IF THIS IS NOT THE CASE, THEN 
C CHANGE THE NEXT STATEMENT APPROPRIATELY. 
C 

IP=6 

WRITE (IP, 998) 

WRITE (IP, 170) 

WRITE(IP,999) 

WRITE (IP, 180) 

WRITE (IP, 999) 

WRITE (IP, 190 )AN1,AN2 

WRITE ( IP , 200 ) XBAR1 , XBAR2 

WRITE ( IP , 210 ) SSQDX1 , SSQDX2 

WRITE ( IP , 220 ) ALP1 , ALP2 

WRITE ( IP , 230 ) BET1 , BET2 

WRITE ( IP , 240 ) RSD1 , RSD2 

WRITE (IP, 999) 

WRITE(IP,250)P 

WRITE (IP, 999) 

WRITE ( IP , 260 ) HPP , NUEQ , GEQ 

WRITE ( IP , 270 ) HPP , NUNEQ , GNEQ 
C 

C COMPUTE THE RATIO OF THE LARGER ESTIMATED RESIDUAL VARIANCE TO 
C THE SMALLER ESTIMATED RESIDUAL VARIANCE 
C 

RATIO= ( RSD1 *RSD1 ) / ( RSD2*RSD2 ) 

AM=ANl-2 . ODO 

AN=AN2-2 . ODO 

IF ( RATIO . GT . 1 . ODO ) G0T0110 

RATI0=1.0D0/RATI0 

HOLD=AM 

AM=AN 

AN=HOLD 
110 WRITE (IP, 999) 

Ml=AM+0 . 1 

Nl=AN+0 . 1 

WRITE ( IP, 280) RATIO 

WRITE (IP, 310) Ml, Nl 

WRITE (IP, 999) 

WRITE (IP, 290 )AML 

WRITE (IP, 999) 

WRITE (IP, 300 )S 

WRITE (IP, 999) 

WRITE (IP, 999) 
C 

C COMPUTE THE COEFFICIENTS OF THE QUADRATIC EQUATION FOR THE EQUAL 
C VARIANCES C/SE 
C 



WRITE (IP, 340) 
WRITE (IP, 999) 
A1=DELBET**2 
A2=S*S*GEQ*GEQ 
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A3= ( 1 . 0D0/SSQDX1 ) + ( 1 . 0D0/SSQDX2 ) 

Bl=2 . ODO*DELALP*DELBET 

B2=2.0D0*A2 

B3= ( XBAR1/SSQDX1 ) + ( XBAR2/SSQDX2 ) 

C1=DELALP**2 

C2=A2 

C3= ( 1 . ODO/AN1 ) + ( 1 . 0D0/AN2 ) + ( XBAR1 *XBAR1/SSQDX1 ) + ( XBAR2*XBAR2/SSQDX 
12) 

A=A1-A2*A3 

B=B1+B2*B3 

C=C1-C2*C3 
C 

C CHECK FOR ZERO COEFFICIENT OF THE QUADRATIC TERM FOR THE EQUAL 
C VARIANCES CASE 
C 

IF(A.NE.0.0D0)G0T0120 

AMID=-C/B 

WRITE(IP,320) 

WRITE (IP, 321) AMID 

G0T0140 
C 

C COMPUTE THE ROOTS OF THE QUADRATIC, AND THE MIDPOINT OF THE ROOTS 
C FOR THE EQUAL VARIANCES CASE 
C 

120 DISC=B*B-4.0D0*A*C 
C 

C CHECK FOR A NEGATIVE DISCRIMINANT FOR THE EQUAL VARIANCES CASE 
C 

IF (DISC . GE . . ODO ) G0T0130 

WRITE (IP, 330) DISC 

G0T0140 
130 R00T1=(-B+DSQRT(DISC) )/(2.0D0*A) 

R00T2= ( -B-DSQRT ( DISC ) ) / ( 2 . ODO*A ) 

AMID=-B/(2.0D0*A) 

WIDTH1=R00T1-R00T2 

WRITE (IP, 350 )R00T2 

WRITE (IP, 360 )R00T1 

WRITE (IP, 370) AMID 

WRITE (IP, 380 )WIDTH1 

WRITE (IP, 999) 

WRITE (IP, 999) 
C 

C COMPUTE THE COEFFICIENTS OF THE QUADRATIC EQUATION FOR THE UNEQUAL 
C VARIANCES CASE 
C 
140 WRITE (IP, 420) 

WRITE (IP, 999) 

A1=DELBET**2 

A2=GNEQ*GNEQ 

A3= ( RSD1 *RSD1/SSQDX1 ) + ( RSD2*RSD2/SSQDX2 ) 

Bl=2 . ODO*DELALP*DELBET 

B2=2 . 0D0*A2 

B3=(RSD1*RSD1*XBAR1/SSQDX1)+(RSD2*RSD2*XBAR2/SSQDX2) 

C1=DELALP**2 

C2=A2 

C3= ( RSD1 *RSD1/AN1 ) + ( RSD2*RSD2/AN2 ) + ( RSD1 *RSD1 *XBAR1 *XBAR1/SSQDX1 ) + 
1 ( RSD2*RSD2*XBAR2*XBAR2/SSQDX2 ) 

A=A1-A2*A3 

B=B1+B2*B3 

C=C1-C2*C3 
C CHECK FOR ZERO COEFFICIENT OF THE QUADRATIC TERM FOR THE UNEQUAL 
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c 
c 
c 
c 

c 
c 
c 



VARIANCES CASE 

IF(A.NE.0.0D0)G0T0150 

AMID=-C/B 

WRITE (IP, 400) 

WRITE (IP, 321) AMID 

RETURN 

COMPUTE THE ROOTS OF THE QUADRATIC, AND THE MIDPOINT OF THE ROOTS 
FOR THE UNEQUAL VARIANCES CASE 

150 DISC=B*B-4.0D0*A*C 

CHECK FOR A NEGATIVE DISCRIMINANT FOR THE UNEQUAL VARIANCES CASE 

IF ( DISC . GE . . ODO ) G0T0160 
WRITE (IP, 410) DISC 
RETURN 
160 R00T1= ( -B4-DSQRT ( DISC ) ) / ( 2 . 0D0*A ) 
R00T2= ( -B-DSQRT ( DISC ) ) / ( 2 . 0D0*A ) 
AMID=-B/(2.0D0*A) 
WIDTH2=R00T1-R00T2 
RATWID=WIDTH2/WIDTH1 
WRITE (IP, 350 )R00T2 
WRITE (IP, 360 )R00T1 
WRITE (IP, 370) AMID 
WRITE (IP, 380 )WIDTH2 
WRITE (IP, 999) 
WRITE (IP, 999) 
WRITE (IP, 390 )RATWID 
RETURN 



C 
C 



170 FORMAT (1H 

1 LINES) 
180 FORMAT (1H 
190 FORMAT (1H 
200 FORMAT (1H 
210 FORMAT (1H 
220 FORMAT ( 1H 
230 FORMAT (1H 
240 FORMAT ( 1H 



,58HSTATISTICAL ANALYSIS OF INTERSECTION OF 2 REGRESSION 



,10HINPUT DATA) 

,9HAN1 = ,D15.8,10X,9HAN2 = ,D15.8) 
,9HXBAR1 = ,D15.8,10X,9HXBAR2 = ,D15.8) 
,9HSSQDX1 = ,D15.8,10X,9HSSQDX2 = ,D15.8) 
,9HALP1 = ,D15.8,10X,9HALP2 = ,D15.8) 
,9HBET1 = ,D15.8,10X,9HBET2 = ,D15.8) 
,9HRSD1 = ,D15.8,10X,9HRSD2 = ,D15.8) 

250 FORMAT (1H ,52H2-SIDED CONFIDENCE INTERVALS WITH PROBABILITY WIDTH 
1,D12.5,13H ARE PROVIDED) 

260 FORMAT ( 1H ,4HTHE ,D12.5,42H PERCENT POINT OF THE T DISTRIBUTION WI 
1TH ,I6,50H DEGREES OF FREEDOM (THE EQUAL VARIANCES CASE) IS ,D12.5 

1) 
270 FORMAT (1H ,4HTHE ,D12.5,42H PERCENT POINT OF THE T DISTRIBUTION WI 

1TH ,I6,52H DEGREES OF FREEDOM (THE UNEQUAL VARIANCES CASE) IS ,D12 

1.5) 
280 FORMAT ( 1H ,98HTHE RATIO OF THE LARGER ESTIMATED RESIDUAL VARIANCE 

1T0 THE SMALLER ESTIMATED RESIDUAL VARIANCE IS ,D15.8) 
290 FORMAT (1H , 35HTHE MAXIMUM LIKELIHOOD ESTIMATOR = ,D15.8) 
300 FORMAT (1H ,4HS = ,D15.8) 
310 FORMAT ( 1H ,79HC0MPARE THE ABOVE RATIO VALUE TO THE PERCENT POINTS 

10F THE F DISTRIBUTION WITH ,I6,5H AND ,I6,19H DEGREES OF FREEDOM) 

320 FORMAT (1H ,86H*****ERR0R — THE COEFFICIENT OF THE QUADRATIC TERM IS 
1 ZERO FOR THE EQUAL VARIANCES CASE) 

321 FORMAT (1H , 57H THE CONFIDENCE INTERVAL HAS ZERO WIDTH A 
1B0UT ,D15.8) 

330 FORMAT (1H ,97H*****ERR0R— NEGATIVE DISCRIMINANT IN THE EQUAL VARIA 
1NCES CASE. THE VALUE OF THE DISCRIMINANT IS ,D15.8) 
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340 FORMAT 
350 FORMAT 
360 FORMAT 
370 FORMAT 
380 FORMAT 
390 FORMAT 

1) 
400 FORMAT 

1 ZERO 
410 FORMAT 

1IANCES 
420 FORMAT 
998 FORMAT 



(1H ,20HEQUAL VARIANCES CASE) 

(1H ,37HL0WER LIMIT OF CONFIDENCE INTERVAL = ,D15.8) 

(1H ,37HUPPER LIMIT OF CONFIDENCE INTERVAL = ,D15.8) 

(1H ,37HMIDP0INT OF CONFIDENCE INTERVAL = ,D15.8) 

(1H ,37HWIDTH OF CONFIDENCE INTERVAL = ,D15.8) 

(1H ,46HTHE RATIO OF THE CONFIDENCE INTERVAL WIDTHS 



,D15.8 



1H ,88H*****ERR0R — THE COEFFICIENT OF THE QUADRATIC TERM IS 
FOR THE UNEQUAL VARIANCES CASE) 
1H ,99H*****ERR0R — NEGATIVE DISCRIMINANT IN THE UNEQUAL VAR 
CASE. THE VALUE OF THE DISCRIMINANT IS ,D15.8) 
1H,22HUNEQUAL VARIANCES CASE) 
1H1) 



999 FORMAT (1H0) 
END 



STATISTICAL ANALYSIS OF INTERSECTION OF 2 REGRESSION LINES 



INPUT DATA 








AN1 


.14000000+002 


AN2 


.11000000+002 


XBAR1 = 


.67500000+002 


XBAR2 = 


-.50000000+001 


SSQDX1 = 


.56875000+004 


SSQDX2 = 


.27500000+004 


ALP1 


.82368597+000 


ALP2 


.83584345+000 


BET1 


.62355168-003 


BET2 


. 23469094-003 


RSD1 


.11459900-003 


RSD2 


.49539027-004 



2-SIDED CONFIDENCE INTERVALS WITH PROBABILITY WIDTH 



.95000+000 ARE PROVIDED 



THE .97500+002 PERCENT POINT OF THE T DISTRIBUTION WITH 
FREEDOM (THE EQUAL VARIANCES CASE) IS .20796+001 

THE .97500+002 PERCENT POINT OF THE T DISTRIBUTION WITH 
FREEDOM (THE UNEQUAL VARIANCES CASE) IS .21098+001 



21 DEGREES OF 
17 DEGREES OF 



THE RATIO OF THE LARGER ESTIMATED RESIDUAL VARIANCE TO THE SMALLER ESTIMATED 

RESIDUAL VARIANCE IS .5351393+001 
COMPARE THE ABOVE RATIO VALUE TO THE PERCENT POINTS OF THE F DISTRIBUTION 

WITH .12000+002 AND .90000+001 DEGREES OF FREEDOM 



THE MAXIMUM LIKELIHOOD ESTIMATOR = .31264354+002 



S = .92500246-004 
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479-206 O - 73 - 7 



EQUAL VARIANCES CASE 



LOWER LIMIT OF CONFIDENCE INTERVAL 
UPPER LIMIT OF CONFIDENCE INTERVAL 
MIDPOINT OF CONFIDENCE INTERVAL 
WIDTH OF CONFIDENCE INTERVAL 



.30804202+002 
.31727843+002 
.31266023+002 
.92364138+000 



UNEQUAL VARIANCES CASE 



LOWER LIMIT OF CONFIDENCE INTERVAL 
UPPER LIMIT OF CONFIDENCE INTERVAL 
MIDPOINT OF CONFIDENCE INTERVAL 
WIDTH OF CONFIDENCE INTERVAL 



. 30865357+002 
.31660331+002 
.31262844+002 
.79497346+000 



THE RATIO OF THE CONFIDENCE INTERVAL WIDTHS 



. 86069493+000 
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